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Abstract 

The VERITAS array of Cherenkov telescopes has carried out a deep observational program 
on the nearby dwarf spheroidal galaxy Segue 1. We report on the results of nearly 48 hours of 
good quality selected data, taken between January 2010 and May 2011. No significant 7-ray 
emission is detected at the nominal position of Segue 1, and upper limits on the integrated flux are 
derived. According to recent studies, Segue 1 is the most dark matter-dominated dwarf spheroidal 
galaxy currently known. We derive stringent bounds on various annihilating and decaying dark 
matter particle models. The upper limits on the velocity-weighted annihilation cross-section are 

(av)95%CL< 10-23 

cm'^ s ^ , improving our limits from previous observations of dwarf spheroidal 
galaxies by at least a factor of two for dark matter particle masses m^ ^ 300 GeV. The lower limits 
on the decay lifetime are at the level of r^^'^*-"^ > 10^^ s. Finally, we address the interpretation 
of the cosmic ray lepton anomalies measured by ATIC and PAMELA in terms of dark matter 
annihilation, and show that the VERITAS observations of Segue 1 disfavor such a scenario. 

PACS numbers: 95.85.Pw, 98.52.Wz, 98.56.Wm, 95.35.+d 
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I. INTRODUCTION 



The compelling evidence for the presence of non-baryonic dark matter in various struc- 
tures in the Universe [l| has motivated numerous efforts to search for dark matter. Among 
many theoretical candidates for dark matter (see iH] for a review of candidates, and experi- 
mental searches), Weakly Interacting Massive Particles (WIMPs) are the most popular and 
well-motivated. A massive thermal relic of the early universe, with a weak scale interaction, 
naturally gives the measured present-day cold dark matter density ficoivih^ = 0.1109 ± 0.0056 

. Candidates for WIMP dark matter are present inmany extensions of the standard model 
of particle physics, such as supersymmetry (SUSY) [3] or theories with extra dimensions j4|. 
In such models, the WIMPs either decay or self-annihilate into standard model particles, 
which can produce a continuum of 7-rays with energies up to the dark matter particle mass, 
or monoenergetic 7-ray lines. Constraints from particle collider experiments are highly 



model-dependent, but general 
of GeV to a few tens of TeV 



y place the mass of such particles in the range of a few tens 
l|]. Indirect searches for dark matter with high energy (HE, 
1 GeV < E < 100 GeV) or very high energy (VHE, 100 GeV < E < 30 TeV) 7-rays thus pro- 
vide a very promising way to test the nature of dark matter. Unlike cosmic ray charged 
particles, HE and VHE 7-rays are free of any propagation effects over short distances (< 
IMpc), and would therefore easily characterize a dark matter source location, spectrum 
and morphology. Such searches are generally conducted using pointed astrophysical obser- 
vations of nearby dark matter overdensities, because the annihilation/decay rate stro ngly 
depends on the dark matter density. Pop ular targets include the Galactic Center js-lll, 



satellite galaxies of the Milky Way |l2l-ll8| , globular clusters [l9|, |20| and clusters of galaxy 



21 



26| . Non-targeted searches are also currently under consideration. They include blind 



searches for dark matter substructures in the galactic halo 
of the galactic and extra-galactic 7-ray diffuse emission 



30 



27H29| and the measurement 
33] . Compared to the targeted 



searches, non-targeted searches are much less affected by the uncertainties in the dark mat- 
ter distribution modeling. However, they can suffer from large uncertainties in the modeling 
of the different backgrounds, especially in the HE 7-ray regime. 

Beyond this well-established picture, additional effects might play an important role in the 
phenomenology of dark matter and the prospects for its detection. Motivated by the recent 



cosmic ray lepton spectra measured by the ATIC [3J], PAMELA 
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36|, H.E.S.S. (37|,l38| 



and Fermi-LAT 



39 



40| experiments in the 1 GeV - 1 TeV energy range, various particle 



physics and astrophysics effects have been suggested which could boost the dark matter 
signal with respect to standard expectations. For example, such particle physics effects in- 



clude the Sommerfeld enhancement to the WIMP annihilation cross-section in t 
matter particle velocity regime 



le low dark 
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44| , or the internal bremsstrahlung effect 



45 



46|, which 



can provide a considerable enhancement to the 7-ray signal at the endpoint of the spectrum. 
The presence of gravitationally-bound substructures within smooth dark matter halos can 
also have a significant impact on the annihilation/decay rate of dark matter particles 
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51| . Such astrophysical enhancements, combined with particle physics en 



lancements, have 



been proposed as an explanation for the cosmic ray lepton anomalies 
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53|. Finally, 



decaying dark matter models have also been suggested to explain the lepton excesses and 
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are good alternatives to annihilating dark matter models 
The dwarf spheroidal galaxies (dSphs) of the Local Group best meet the criteria for a clear 
and unambiguous detection of dark matter. They are gravitationally-bound objects and 
are believed to contain up to (9(10^) times more mass in dark matter than in visible mat- 
;er , making them widely discussed as potential targets for indirect dark matter detection 
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-|64|. Dwarf spheroidal galaxies are believed to be the remnants of dark matter 
halos, which contributed to the formation of Milky Way-sized galaxies during hierarchi- 
cal clustering in structure formation scenarios. As opposed to the Galactic Center, and 
possibly globular clusters 65|, l66|], they are environments with a favorably low astrophys- 
ical 7-ray background. Neither astrophysical 7-ray sources (supernova remnants, pulsar 
wind nebulae, etc) nor gas acting as target material for cosmic rays have been observed 
in these systems 67|, |68|]. Furthermore, their relative proximity and high galactic latitude 
make them the best astrophysical targets for high a signal-to-noise detection. With star 
velocity dispersions of the order of 10 kms"^ [69], dSphs are ideal laboratories for testing 
a possible velocity-dependence of the dark matter annihilation cross-section (for instance 
the Sommerfeld enhancement predicts (crv) oc 1/v). Like the Milky Way halo, dSphs are 
thought to harbor a population of substructures that could possibly boost their overall dark 



matter luminosity. However, recent simulations and analytic calculations show that 
these objects, the expected astrophysical boost is less than a factor of ten 
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63|. 



or 



The sensitivity improvement of the latest infrared/optical sky surveys has doubled the 
known number of Milky Way satellites in the past few years (for instance, see the Sloan 
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Digital Sky Survey (SDSS) recent discoveries [70|]), renewing the interest in these objects 
as promising targets for indirect dark matter searches. Segue 1 is one of these new Milky 
Way satellites, an ultra-faint dSph discovered in 2006 as an overdensity of resolved stars 
in the SDSS j70|. It is located at a distance of 23 ± 2 kpc from the Sun at (RA,Dec) 
= (10^^07^03. 2^, 16°04'25"), well above the galactic plane. Because of its proximity to the 
Sagittarius stream, the nature of the Segue 1 overdensity has recently been disputed, with 
some authors arguing that it was a tidally disrupted star cluster originally associated with 
the Sagittarius dSph [7l|. However, a kinematic study of a larger member-star sample (66 
stars compared to the previous 24-star sample) has recently confirmed that Segue 1 is an 
ultra- faint Milky Way satellite galaxy 72|. According to a study of its star kinematics. 
Segue 1 is probably one of the most dark matter-dominated dSph and is often highlighted 
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72 



73| . Segue 



75| in its survey 



as the most promising dSph target for indirect dark matter searches 
1 has been observed in the HE 7-ray regime by the Fermi-LAT satellite 
observation mode. Although no data analysis has been published yet by the Fermi-LAT 
collaboration, dedicated searches for dark matter using the first 9 months of public Fermi- 



LAT data on Segue 1 have been carried out by several authors 18 
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No 7-ray signal 



was discovered in any of these analyses. The resulting upper limits on the dark matter 



velocity-weighted annihilation cross-section are at the level of (av) ~ 10 



-21 



10 



-23 



cm s 



in the 10 GeV - 1 TeV WIMP mass range. In the VHE band, the MAGIC collaboration has 
recently conducted a search for dark matter annihilation in Segue 1, analyzing a 30-hour 
dataset taken in single telescope mode 76[. No VHE 7-ray signal was discovered, giving 
upper limits on (av) at the level of ~ 10^22 _ ^Q-23 ^^3 ^-i ^j^^ ^qq Q^y _ 2 XeV WIMP 
mass range. 

This paper reports on extensive observations of Segue 1 conducted by the Very Energetic 
Radiation Imaging Telescope Array System (VERITAS). After describing the VERITAS 
instrument, the observations and the data analysis in section [Tll we extract integrated flux 
upper limits assuming various types of spectra in section IIIII In section IIVI bounds on 
annihilating and decaying dark matter are derived. Section |V] addresses the interpretation 
of the recent cosmic ray lepton anomalies in terms of dark matter annihilation and presents 
the VERITAS 7-ray constraints using the Segue 1 data. Finally, section |VT] is devoted to 
our conclusions. 
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II. OBSERVATIONS AND ANALYSIS 



VERITAS is an array of four 12-meter imaging atmospheric Cherenkov telescopes 
(lACTs) located at the base camp of the F. L. Whipple Observatory in southern Ari- 
zona (31°.68 N, 110°.95 W, 1.3 km above sea level). Each VERITAS telescope consists of 
a large optical reflector which focuses the Cherenkov light emitted by particle air showers 
onto a camera of 499 photomultiplier tubes. The array has been fully operational since 
September 2007. The large effective area (~ 10^ m^), in conjunction with the stereoscopic 
imaging of air showers, enables VERITAS to be sensitive over a wide range of energies 
(from 100 GeV to 30 TeV) with an energy and angular resolution of 15%-20% and 0.1° per 
reconstructed 7-ray, respectively. VERITAS is able to detect a point source with 1% of the 
Crab Nebula flux at a statistical significance of 5 standard deviations above background 
(5o") in approximately 25 hours of observations. For further details about VERITAS, see, 
e.g., H- 

Observations of the Segue 1 dSph were performed between January 2010 and May 2011. 
The data used for this analysis only includes observations taken under good weather and 
good data acquisition conditions. The total exposure of the dataset, after quality selection 
and dead time correction, amounts to 47.8 hours, and is the largest reported so far for any 
dSph observations conducted by an array of lACTs. The mean zenith angle of the observa- 
tions is ~ 20°. The observations were conducted using the wobble pointing strategy, where 
the camera center is offset by 0.5° from the target position. The wobble mode allows for 
simultaneous background estimation and source observation, thus reducing the systematic 



uncertainties in the background determination 78 1. 



The data were reduced using standard VERITAS calibration and analysis tools 79|. After 



calibration of the photomultiplier tube gains 



80| , images recorded by each of the VERI' 



telescopes are characterized by a second-moment analysis giving the Hillas parameters 



AS 



A stereoscopic analysis combining each telescope's image parameters is then used to recon- 



struct the 7-ray arrival direction and energy 82| . We applied selection criteria (cuts) on the 
mean-reduced-scaled length and mean-reduced-scaled width parameters (see 83| for a full 
description of these parameters) to reduce the hadronic cosmic ray background. The cuts 
for 7-ray/hadron separation were optimized a priori for a source with a 5% Crab Nebula- 
like flux. Additionally, an event is accepted as a 7-ray candidate if the integrated charge 



recorded in at least two telescopes is > 90 photoelectrons, which effectively sets the analysis 
energy threshold to 170 GeV. Finally, a cut on 9, the angle between the target position 
and the reconstructed arrival direction, is applied to the 7-ray candidates and defines the 
signal search region (0^ < 0.015 deg^ in our analysis). After 7-ray selection, the residual 
background was estimated using the ring background technique |i8^]. The ring background 
method computes the background for each position in the field of view using the background 
rate contained in a ring around that position. Two circular regions, of radius 0.2° centered 
on the target position and of radius 0.3° centered on the bright star ?7-Leonis (with apparent 
magnitude in the visible band My = 3.5, and located 0.68° from the position of Segue 1), 
were excluded for the background determination. 

The analysis of the data resulted in the selection of Nqn = 1082 7-ray candidates in the 
signal search region and Nqff = 12479 background events in the background ring region, 
with a normalization factor a = 0.084, resulting in 30.4 excess events. The corresponding 
significance, calculated according to the method of Li & Ma 85|, is 0.9 cr. No significant 
7-ray excess is found at the nominal position of Segue 1, nor in the whole field of view, 
as shown by the significance map on Figure [H The large depletion area, with negative 
significances, corresponds to the bright star //-Leonis. 



<u 
o 

4 n 
o 

'E 

u> 

2 ■« 

-2 
-4 
"6 

right ascension [hours] 

^ J2000 

FIG. 1. Significance map obtained from the VERITAS observations of Segue 1 after 7-ray selection 
and background subtraction. The black cross indicates the position of Segue 1. The black circles 
correspond to the two exclusion regions used for the background determination. See text for further 
details. 
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Given the absence of signal, one can derive upper limits (ULs) on the number of 7-rays in 
the source region. The computation of statistical ULs can be done following several meth- 
ods, each relying on different assumptions. The bounded profile likelihood ratio statistic 



developed by Rolke et al. [86| is used in our analysis. As discussed in the following sections, 
these ULs will serve for the computation of integrated flux ULs and for constraining some 
dark matter models. To make the computation of integrated flux ULs robust, we define a 
minimum energy, above which the energy reconstruction bias is less than 5%. The energy 
reconstruction bias as a function of the reconstructed 7-ray energy has been studied with 
Monte Carlo simulations. The 7-ray selection cuts used in this analysis set the minimum 
reconstructed energy to Emin = 300 GeV. The ULs on the number of 7-rays computed with 
the Rolke prescription are displayed in Table HI along with the analysis results. 



Live time (min) 


Emin (GeV) 


N;-(E > E^i„) 


Significance 


N95%CL(E>E^i,) 


2866 




30.4 


0.9 


135.9 


2866 


300 


31.2 


1.4 


102.5 



TABLE L Analysis results of the VERITAS observations of Segue 1. N^'''=(E > Emin) is the number 
of excess events in the signal search region with energies E > Emin, after background subtraction. 

N95%CL(e>E 

min) is the 95% confidence level (CL) upper limit on the numb er of 7-rays with 



energies E > Emin in the signal search region, computed according to the Rolke 



86( prescription. 



III. FLUX UPPER LIMITS 

The analysis of the data did not show any significant excess over the background at the 
nominal position of Segue 1. The ULs on the number of 7-rays in the signal search region 
can then be converted to ULs on the integral 7-ray flux. The number of 7-rays detected 
by an array of lACTs above a minimum energy Emin is related to the source integral flux 
$^(E > Emin) by: 

TOO A f E) '^^^ dE 

N,(E > Emin) = Tobs X ' dN. 'g >< "^^(^ > Emin), (l) 

Jtjmin dE 

where Tobs is the observation time, dN^/dE the assumed source differential energy spectrum 
and v4efr(E) is the instrument effective area. The effective area v4eff(E) is the instrument 



response function to the collection of 7-rays of energy E, and it depends on the zenith angle 
of the observations, the offset of the source from the target position and the 7-ray selection 
cuts. In the next two subsections, we consider two assumptions for the differential 7-ray 
spectrum: the case of a generic power-law spectrum, which describes well the TeV energy 
spectra of standard astrophysical sources and the case of a 7-ray spectrum resulting either 
from the annihilation or the decay of WIMP dark matter. 



A. Upper limits with power-law spectra 

Table Hn shows the integral flux ULs above Emin 
law spectrum: 

dN. 

oc E 



dE 



300 GeV for the assumption of a power- 



(2) 



where F is the spectral index. The spectral indices have been varied over the range 
r = 1.8 — 3.0. The ULs on the integrated flux do not depend on the flux normalization 
(see eq. [1]), but they do depend on the spectral index. A harder power-law spectrum 
provides a more constraining upper limit. The ULs on the integral flux reported in table ITTl 
are at the level of 0.5% of the Crab Nebula integral flux. 



Spectral index 

r 


^95%CL(E > 300 GeV) 
[10-^2 cm^^s-i] 


1.8 


7.6 


2.2 


7.7 


2.6 


8.0 


3.0 


8.2 



TABLE II. 95% CL ULs on the integrated 7-ray flux above Emin = 300 GeV from the VERITAS 
observations of Segue 1, for power-law spectra with various spectral indices. For comparison, 1% 
of the integrated Crab Nebula flux above Emin = 300 GeV is 1.5 x 10~^^ cm~^ s~^. 
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B. Upper limits with Dark Matter 7-ray spectra 

Since Segue 1 is dark matter-dominated, one can derive ULs on the integrated flux as- 
suming that the dominant source of 7-rays is dark matter annihilation or decay. The 7-ray 
differential energy spectrum from dark matter particle annihilation or decay depends on the 
dark matter model, and especially on the branching ratios to the final state particles. In 
almost every channel (excepting the e"^e~ and fi~^fi~ channels), the 7-ray emission mostly 
originates from the hadronization of the final state particles, with the subsequent produc- 
tion and decay of neutral pions. Annihilation/decay to three different final state products 
is considered independently of the dark matter model, in each case with a 100% branching 
ratio: W"^W~, bb and r+r^. For each channel, the 7-ray spectrum has been simulated 
with the particle physics event generator PYTHIA 8.1 [87]. As shown by the left panel 
of Figure [2], the bb and r+r~ channels encompass a wide range of dark matter 7-ray an- 
nihilation/decay spectra and give an idea of the uncertainties related to the dark matter 
particle physics model. The Z^Z*^ channel gives a 7-ray spectrum very similar to the W"'"W~ 
channel and is not considered here. The right panel of Figure [2] shows the 95% CL ULs on 
the integrated flux above Emin = 300 GeV as a function of the dark matter particle mass. 
For m^ < 800 GeV, the most constraining ULs are obtained for the bb channel because the 
bb spectrum features a small "bump" at high x = E/m^ values (see Figure [2] left), which 
makes it the hardest among all the considered spectra. Increasing the dark matter particle 
mass, the x lower limit over which the 7-ray spectrum is integrated extends down to lower 
values, making the bb 7-ray spectrum on average softer (and, inversely, the r^r^ spectrum 
harder). Above a dark matter particle mass of 800 GeV, the r+r^ spectrum then gives the 
most constraining integrated flux ULs. The integrated flux ULs range between 0.3% and 
0.7% of the Crab Nebula integral flux, depending on the dark matter particle mass. 



IV. DARK MATTER BOUNDS 

The absence of signal at the position of Segue 1 can be used to derive constraints on 
various dark matter models. Two different scenarios, in which the modeling of the 7-ray 
flux slightly differs, are considered: the case of annihilating dark matter and the case of 
decaying dark matter. 
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FIG. 2. Left: dark matter annihilation/decay spectra for four different final state products 



871]. The spectra are plotted in 



(W+W-, Z^Z", bb and r+r"), extracted from PYTHIA 8.1 
the dN/dx representation, where x = E/m^ (or x = 2E/m^ for decay spectra). Right: 95% CL 
ULs on the integrated 7-ray flux above Emin = 300 GeV from the VERITAS observations of Segue 
1 considering dark matter particle annihilation/decay for three different channels: W^W^, bb 
and r+r^. For comparison, 0.5% of the integrated Crab Nebula flux above Emin = 300 GeV is 
7.5 X lO^^^cm-^g-i. 



A. 7-rays from dark matter annihilation or decay 



The differential 7-ray flux from the annihilation of dark matter particles of mass m^, in 
a spherical dark matter halo is given by a particle physics term multiplied by an astrophysical 



term 



59|: 



d<l>^ 



(AO,E) 



1 (av) dN^ 



4:71 2 ml dE 

A. 



(3) 



The particle physics term contains all the information about the dark matter model: the 
mass of the dark matter particle m^, the 7-ray differential energy spectrum from all final 
states weighted by their corresponding branching ratios dN^/dE, and its total velocity- 
weighted annihilation cross-section (av). The astrophysical factor J(Ai7) (sometimes called 
the dark matter annihilation luminosity) is the square of the dark matter density integrated 
along the line of sight, s, and over the solid angle AQ 59| : 



(4) 
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where the upper and lower hne of sight integration bounds depend on the distance d and 
the tidal radius rt of the target: 



Smax.min = dcos(6') ± ^r? - d^ siu [6). (5) 

The solid angle AQ = 27r (1 — cos(0max)) is given here by the size of the signal search region 
defined previously in our analysis, i.e. 9^^^ = 0.015 deg^. The estimate of the astrophysical 
factor requires a model of the Segue 1 dark matter profile. Motivated by results from N-body 



simulations, an Einasto profile 88N90| is used: 

p,(r)=p.e--[(^/-)^^"-l, (6) 
with the scale density, the scale radius and the index n respectively being ps = 1.1 x 10^ Mq kpc" 



rs = 0.15kpc and n = 3.3 9l|]. The value of the tidal radius changes the astrophysi- 
cal factor by less than 10%. We adopt a value of 500 pc, which is the median trun- 
cation radius of the Via Lactea II simulation subhalos presenting similar characteristics 



to the Segue 1 dark matter halo [49|, |72|. Having these parameters in hand, the value 
of the astrophysical factor within the solid angle subtended by our on-source region is 
J(Ar2) = 7.7 X 10^*^ GeV^^ cm~^ sr. The systematic uncertainties on the astrophysical fac- 
tor resulting from the fit of the Segue 1 dark matter distribution to an Einasto profile are 



73|. 



less than an order of magnitude at the la level 
In the scenario where the dark matter is a decaying particle, the expression of the differential 
7-ray flux slightly differs and can be obtained with the following substitutions: 

• in eq. [3] ((Tv)/2m^ — )• F/m^, where F = r^^ is the inverse of the dark matter particle 
decay lifetime. 

• in eq. IH p^- 

Because particle decay is a one-body process, the astrophysical factor now depends on the 
dark matter density only. Adopting the same profile parameters, the value of the astrophys- 
ical factor in decaying dark matter scenarios is J(Ai7) = 2.6 x 10^^ GeV cm~^ sr. As in the 
annihilating dark matter case, the uncertainties on the astrophysical factor computed for an 
Einasto dark matter profile are less than one order of magnitude at the la level [tJ]. 
The total number of 7-rays above a minimum energy Emin expected in the signal source re- 
gion is the integral of the differential 7-ray flux (eq. [3]), taking into account the instrument 
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response to the collection of 7-rays: 

N^(E > E^i„) = Tobs X / Aff(E)--?^ dE (7) 

This equation is equivalent to eq. [T] used for the calculation of integral flux upper limits 
(see section IIIII) if one assumes that the only source of 7-rays in Segue 1 is dark matter 
annihilation or decay. From eq. [7] and taking Emin = 0, the ULs on the number of 7- 
rays previously derived in section |IT1 (see table [T]) can either be translated into ULs on the 
velocity-weighted annihilation cross-section (crv) or lower limits (LLs) on the decay lifetime 
r. 



B. Upper limits on the annihilation cross-section 

In this section, ULs on the dark matter particle annihilation cross-section (crv) are derived 
independently of any dark matter models, considering five different final states with 100% 
branching ratios: W^W~, bb, r"''r~, e"'"e~ and . The W"''W^, bb and r+r~ differential 



7-ray spectra have been simulated with PYTHIA 8.1 [87| (see section [Tll B|) and are displayed 
on the left panel of Figure O The final state radiation (FSR) e"'"e~ and yU^/i— channels are 
motivated by the recent anomalies measured in the cosmic ray lepton spectra (see section 
|V|) . The FSR spectrum for lepton channels has been computed analytically in (45| for the 
e"'"e~ case, and is given by: 



where a is the fine structure constant, mf the fermion particle mass, x = E/m^ and 



/3f = y 1 — mf/m^ X (1 — x). This analytical formula has been cross-checked against PYTHIA 
simulations and has given good agreement for both the e"'"e~ and fi'^fi~ channels. In addi- 
tion to the FSR contribution of the muons, the contribution of the radiative muon decay 



92|. 



fi^ — )■ e~ z/^ z/e 7 and yU+ — )■ e"*" Ue 7 has been included in the /U~ 7-i'ay spectra 
Figure [3] shows the 95% CL exclusion curves on (crv) as a function of the dark matter 
particle mass for the five channels considered above, using eq. [3] and eq. [71 For the 
W^W~ channel, the 95% CL UL on the velocity-weighted annihilation cross-section is 
(crv)^^^"*-"^ < 8 X 10~^^cm^s~^ at 1 TeV. This limit is the most constraining reported so 
far for any dSph observations in the VHE 7-ray band. The bb and r+r~ exclusion curves 
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illustrate the range of uncertainties on the (av) ULs from the dark matter particle physics 
model. Concerning the lepton channels e"'"e~ and yU'^/i", the limits are at the level of 
10^^"^ cm'^ s~^ at 1 TeV. The current ULs on (av) are two orders of magnitude above the 
predictions for thermally produced WIMP dark matter. 




m-; [GeV] m-, [GeV] 

FIG. 3. 95% CL ULs from the VERITAS observations of Segue 1 on the WIMP velocity-weighted 
annihilation cross-section (cv) as a function of the WIMP mass, considering different final state 
particles. The grey band area represents a range of generic values for the annihilation cross-section 
in the case of thermally produced dark matter. Left: hadronic channels W'^W^, bb and t'^t~ . 
Right: leptonic channels e~''e~ and . 

C. Lower limits on the decay lifetime 

If we assume that dark matter is a decaying particle, LLs on the lifetime of dark matter 
can be derived. In decaying dark matter scenarios, the dark matter particle can either 
be bosonic or fermionic. The LLs are computed using eq. [7] and making the appropriate 
substitutions to eq. [3], as explained in section IIV A[ For bosonic dark matter particles, the 
same channels as in the annihilating dark matter case are considered: W"''W^, bb, r+r^, 
e"'"e~ and fi'^fi". The decay spectra are the same as those used for the annihilating dark 
matter bounds (see right panel of Figure [21 and eq. [H]), making the substitution for the 
scaled variable x — > 2x, or equivalently m^ — )■ m^/2. The left panel of Figure H] shows the 
95% LLs on the decay lifetime r for the five channels mentioned above. The limits peak at 
the level of r ~ 10^^ — 10^^ s, depending on the dark matter particle mass. 
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In the case of fermionic dark matter, the dark matter particle decays to different channels. 
The following channels are considered for the exclusion limits: W^£^ (where i = e, n, r) and 
TPv. Each corresponding 7-ray spectrum is a combination of the same spectra used in section 
IIVBI for the computation of limits on annihilating dark matter. They have been simulated 
with the PYTHIA 8.1 package |^]. The right side of Figure S] shows the corresponding LLs 
on r, which peak in the range r ~ 10^*^ — 10^^ s. 

Interestingly, decaying dark matter models have been suggested to be good alternatives to 
annihilating dark matter models for explaining the ATIC 3J] and PAMELA [35] lepton 
anomalies, because they circumvent model-building issues such as the ad hoc addition of 
boost factors (e.g Sommerfeld enhancement, and/or astrophysical boost factors, see section 
|V|) . Spectral fits to the Fermi-LAT and PAMELA data prefer channels involving hard lepton 
spectra such as /i^yU^, r+r~ and W^/i^, with a dark matter particle in the mass range 2-5 



TeV and with a lifetime of the order of r ~ 10^^ s 5J, |57|, |58|]. The VERITAS limits on the 
dark matter particle decay lifetime are at least an order of magnitude away from the best 
fits to the Fermi and PAMELA data (see Figure HI) and do not rule out these models. 




10* 

m [GeV] 



. W*e" or W e* 

> W*n" or W'|i* 

> W*x- or W x* 
>Z°v 




10" 
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FIG. 4. 95% CL LLs from the VERITAS observations of Segue 1 on the decay lifetime as a function 
of the dark matter particle mass. Left: Bosonic dark matter decaying to two identical particles: 
W^W^, bb, r+r~, e"'"e~ and The black star and the black cross denote the best fits to 

the Fermi and PAMELA data considering the and r^r^ channels, respectively, and are 



taken from 



58l |. Right: Fermionic dark matter decaying to two different particles: W^^^ (where 
e, n, r) and TPv. The black triangle indicates the best fit to the Fermi and PAMELA data 



considering the channel W^//^, taken from 



58|. 
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V. TESTING THE DARK MATTER INTERPRETATION OF THE COSMIC RAY 
LEPTON ANOMALIES 



The excess in the cosmic ray electron spectrum measured by the ATIC collaboration 



34l ]. and the unexpected rise of the positron fraction 35| in conjunction with the lack of 



anti-proton excess ,94| reported by the PAMELA collaboration, have received considerable 
attention over the past few years. Even if these features can easily be explained by con- 



ventional astrophysical sources 



93 



95| , the dark matter interpretation has been extensively 



studied and has led to numerous dark matter models. In any dark matter annihilation inter- 
pretation, the cosmic ray lepton excesses require a dark matter particle mostly annihilating 
into leptons and with a significant (9(10^ — 10^) boost to the thermal freeze-out annihilation 
cross-section (crv) ~ 3 x lO^^^cm^s^^. The substructures residing in the Milky Way dark 
matter halo do not provide the necessary boost factor [4^ and are unlikely to be responsible 



for the lepton anomalies 



961]. However, models including a Sommerfeld enhancement to the 



annihilation cross-section naturally solve this issue. In this section, we use the VERITAS 
Segue 1 observations to test such models and, more generally, to derive limits on the boost 
factor in a model-independent way. 

A. Models with a Sommerfeld enhancement 



The Sommerfeld enhancement is a non-relativistic quantum effect arising when two dark 
matter particles interact through an attractive potential 43|], mediated by a particle 0. 
The Sommerfeld correction S (or Sommerfeld boost) is velocity-dependent (S ~ 1/v) and 
modifies the product of the annihilation cross-section and the relative velocity: 



av = S(v, m^, m0, a) x (crv)o. 



(9) 



where (ctv)o is the WIMP annihilation cross-section times its relative velocity at thermal 
freeze-out. As shown by eq. [9l the Sommerfeld correction also depends on the dark matter 
particle mass m^, the mass m^ of the particle mediating the attractive potential, and its 
coupling a to the dark matter particle. Depending on the mass and the coupling of the 
exchanged particle, the Sommerfeld enhancement can exhibit a serie of resonances f or sp ecific 
values of the dark matter particle mass, giving very large boost factors up to 10^ 43|]. The 
Sommerfeld enhancement is of particular interest for cold dark matter halos like dSphs, where 
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the mean dark matter velocity dispersion can be as low as a few kms~^. The computation 
of the Sommerfeld enhancement for a relative velocity v, a dark matter particle mass m^, a 
mediator mass m,^ and a coupling a is usually performed through the numerical resolution 
of the Schrodinger equation, modeling the attractive potential with a Yukawa potential 



43(1 . Instead of numerically solving the Schrodinger equation for each set of parameters 



v,m^,m(^,a), we use an analytic solution by approximating the Yukawa potential with the 



Hulthen potential 
numerical solution 



97 
98 



98| . The analytic solution has been shown to closely match the 
. In order to calculate the Sommerfeld enhancement in Segue 1, one 
has to average the Sommerfeld boost factor over the dark matter relative velocity distribution 
f(v): 

S(m^,m<^,a) = y" S(v, m^, m<^, a)f(v) dv. (10) 
The velocity-weighted annihilation cross-section in the dark matter halo is then given by 

(o-v) = (crv)o X S(m^,m^,a). (11) 



Following [73|], the Segue 1 dark matter relative velocity distribution is assumed to be 
Maxwellian, i.e. the dark matter gas is thermalized and at equilibrium, with a mean relative 
velocity dispersion of vq ~ 6.4kms~^. 



eld enhancement to the an- 
assumes that the dark 



In this section, we focus on two models comprising a Sommer: 
nihilation cross-section. The first model (hereafter model I, 
matter particle is a wino-like neutralino x*^, arising in SUSY extensions of the standard 
model. To circumvent the helicity suppression of the annihilation cross-section into light 
leptons, the neutralino can oscillate with charginos ^ which themselves can preferentially 
annihilate into leptons. The transition to a chargino state is mediated by the exchange of 
a Z*^ boson (mgo ~ 90 GeV, a ~ 1/30), leading to a Sommerfeld enhancement. The second 



model (hereafter model II) introduces a new force in the dark sector [4^. The new force 
is carried by a light scalar field predominantly decaying into leptons and with a mass 
(9(1 GeV) and coupling to standard model particles chosen to prevent the overproduction 
of antiprotons. In such models, dark matter annihilates to a pair of scalar particles, with 
an annihilation cross-section boosted by the Sommerfeld enhancement. The coupling a of 
the light scalar particle to the dark matter particle is determined assuming that XX ~^ 00 
is the only channel that regulates the dark matter density before freeze-out {9^. 
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Figure [5] shows the VERITAS constraints for each of these models, derived with the ob- 
servations of Segue 1. The dashed curves show the 95% CL exclusion limits without the 
Sommerfeld correction to the annihilation cross-section, whereas the solid curves are the 
limits to the Sommerfeld enhanced annihilation cross-section. The left panel of Figure 
shows the constraints on model I, for the annihilation of neutralinos into W+W^ through 
the exchange of a Z° boson. The Sommerfeld enhancement exhibits two resonances in the 
considered dark matter particle mass range, for m^ ^ 4.5 TeV and m^ ~ 17TeV, respec- 
tively. VERITAS excludes these resonances, which boost the annihilation cross-section 
far beyond the canonical (av) ~ 3 x 10~^^ cm'^ s~^. The right panel of Figure O shows the 
VERITAS constraints on model II, for a scalar particle with mass m^ = 250 MeV. The 
Sommerfeld enhancement exhibits many more resonances, located at different dark matter 
particle masses and with different amplitudes with respect to model I, because the coupling 
and mass of the exchanged particle differ. Two channels in which the scalar particle decays 
either to e"'"e~ or fi~^^~ have been considered. VERITAS observations start to disfavor 
such models, especially for the e"'"e~e"'"e~ channel where some of the resonances are beyond 
((Tv) ~ 3 X 10^^^ cm'^ s~^. This result holds for particle masses up to a few GeV. 



B. Model-independent constraints on the boost factor 

In the previous section, we have explicitly constrained the Sommerfeld boost factor to 
the annihilation cross-section in the framework of two interesting models. Here, an example 
of model-independent constraints on the overall boost factor Bp (particle physics and/or 
astrophysical boost) as a function of the dark matter particle mass is presented. The con- 
straints are then compared to the recent cosmic ray lepton data. 

Following {99! , we assume that dark matter annihilates exclusively into muons with an anni- 
hilation cross-section (crv) = 3 x 10~^^cm^s~^. In such a case, we use the dashed exclusion 
curve of Figure [3] (right) to compute 95% limits on Bp. Figure E] shows the 95% CL ULs on 
the overall boost factor Bp. The blue and red shaded regions are the 95% CL contours that 
best fit the Fermi-LAT and PAMELA e"'"e~ data, respectively. The grey shaded area shows 



the 95 % CL excluded region derived from the H.E.S.S. e"'"e data 99|. The black dot is an 



example of a model which simultaneously fits well the H.E.S.S., PAMELA and Fermi-LAT 
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FIG. 5. 95% CL exclusion curves from the VERITAS observations of Segue 1 on {av)/S as a 
function of the dark matter particle mass, in the framework of two models with a Sommerfeld 
enhancement. The expected Sommerfeld enhancement S applied to the particular case of Segue 1 
has been computed assuming a Maxwellian dark matter relative velocity distribution. The grey 
band area represents a range of generic values for the annihilation cross-section in the case of 
thermally produced dark matter. Left: model I with wino-like neutralino dark matter annihilating 
to a pair of W^W~ bosons. Right: model II with a 250 MeV scalar particle decaying into either 
e^e~ or n^n^- See text for further details. 



data. The VERITAS VHE 7-ray observations of Segue 1 rule out a significant portion 
of the regions preferred by cosmic ray lepton data. However, the electron and positron 
constraints depend on the cosmic ray propagation model, especially on the electron energy 



loss parameter tq 



99| . The Klein-Nishina suppression of the inve rse Compton loss rate can 



lOOj . In such a case, the Fermi 



significantly alter this parameter at energies above 100 GeV 
and PAMELA 95% CL contours would scale down, and the VERITAS limits would then be 
relaxed. 



VI. CONCLUSIONS 

The ground-based VHE 7-ray observatory VERITAS has started an extensive observa- 
tion campaign toward the nearby Segue 1 dSph, one of the most dark matter-dominated 
satellite galaxies currently known. With nearly 48 hours of exposure, we derived strong flux 
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FIG. 6. 95% CL exclusion limits on the overall boost factor Bp as a function of the dark matter 
particle mass from the VERITAS observations of Segue 1, assuming that the dark matter par- 
ticles annihilate exclusively into fi^fi". The shaded areas are the 2 a contours derived from fits 
to the H.E.S.S., PAMELA and Fermi-Lat data. The black dot is an example of a model that 



simultaneously fits well the H.E.S.S., PAMELA, and Fermi-LAT data. See 



991] for further details. 



ULs and constraints on the annihilation cross-section or decay life time of a hypothetical 
WIMP, independently of any dark matter models. The reported integral flux ULs at the 
95% CL are at the level of 0.5% of the Crab Nebula flux above a minimum energy of 300 
GeV. The corresponding limits on the velocity-weighted annihilation cross-section are in the 
range (av)^^^"'-^^ < 2 — 9 x 10^^^ cm^ s^^, depending on the annihilation channel considered. 
These are the most constraining limits reported so far with any dSph observations conducted 
in the VHE 7-ray band. The limits are complementary to those provided by the Fermi-LAT 



collaboration [17|, but are at least two orders of magnitude away from the canonical value 
(av) ~ 3 X 10~^^cm^s~^. Bounds on the lifetime of decaying dark matter have also been 
derived, with LLs in the range r^^^°^^ > 10^^ — 10^^ s, an order of magnitude below the 
models that best fit the electron and positron excesses recently measured in cosmic ray 
spectra. Finally, the VERITAS Segue 1 data have also been used to test the dark matter 
interpretation of the cosmic ray lepton anomalies. The Segue 1 data disfavors annihilating 



98|. 



dark matter models with a Sommerfeld enhancement, confirming the results of 15| and 
Furthermore, the VERITAS observations start to exclude scenarios where the dark matter 
annihilates preferentially into a /x^/x^ pair, which is the favored scenario for good fits to the 
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H.E.S.S., PAMELA and Fermi-LAT electron and positron data 99|. 

The uncertainties on the dark matter hmits derived throughout this paper mostly come from 
the modeling of the Segue 1 dark matter density profile. As opposed to the classical dSphs, 
the lack of a high statistics star sample for Segue 1 prevents an accurate modeling of the 
dark matter distribution. Assuming an Einasto profile, the systematic uncertainties in the 
dark matter profile modeling can change the astrophysical factor, and hence they scale up 



73| . The systematic uncertainties could 



or down the limits by a factor of 4 at the 1 o level 
even be larger if the Segue 1 dark matter profile is compatible with a cored profile. Future 
spectroscopic surveys might increase the Segue 1 star sample and reduce the uncertainties 
on its dark matter content. 
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